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ABSTRACT 


Properties of two alternative procedures to the Jack- 
knife Point and Confidence Interval Estimation Procedure of 
Gaver and Chu have been studied. They are called the Loq- 
Normal Likelihood Procedure (LNLJ) and the Moment Procedure 
(MP). These two procedures were investigated and compared 
with the Jackknife Point and Confidence Interval Availabil- 
meyeostaMation Procedure. Numerical results from sinula- 


tions are presented in this report. 
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he OVERVIEW 

Availability is the measure of equipment effectiveness 
that relates reliability and maintainability to osperational 
readiness. In some cases, availability and operational 
readiness have been considered the same. In general, these 
are all requirements which must be considered during the 
design and operation of a system, and, thus, must be quanti- 
tatively evaluated. 

There are a variety of — of expressing availability. 
In general, availability relates up time (reliability- 
related) to down time (maintainability- related), and it may 
be defined as the ratio between the time the system is capa- 
ble of performing its mission to the total time the systen 
is in operational demand. Alternatively, it is a measure of 
the probability that a single system is "up" or available at 
a possibly "random" point in time when its services are 
needed. 


The three expressions of availability of greatest con- 


Sern are ~—(1) inherent availability (Ai), (2) operational 





mugatability (Ao), and (3) achieved availability (Aa). See 


Rise and Bjorklund{[ Ref. 1]. In this study, we will be con- 
cerned with only the inherent availability. Often, inherent 
availability, which is a hardware oriented measure, is spe- 
cified and required within a maintainability contract 
reguirement. 

A matter of primary concern is that of estimating avail- 
ability from observations on system up times and down times. 
The estimating process should furnish both a point, or sin- 
gle number, estimate and also somé neasure of the stability 
of the estimate such as a standard error, or confidence lin- 
its.The latter problem is more difficult than ths former. 

In the paper by Gaver and Chuf Ref. 2}, it has been 
demonstrated that the jackknife technique can be useful for 
confidence interval estimation of system availability.¥for 
instance, let 


(0) 
 _———— (1.1) 
2 (U0) +2(D) 


be the long run system availability of a single unit system 
that changes from up to down states in accordance with a 
two-state semi-Markov (or more general) process for which up 


{down} state expected duration is E{U)[5(D) J]. Then it is 





possible +o successfully jackknife the naive point estimator 
that involves replacing expectations by the corresponding 


sample means, 


Qt 


A= ————— Woe 

U + D 

after initial logistic transformation. There are, however, 
particular parametric families of distributions, popular as 
models for summarizing up and down time data, for which the 
most statistically efficient estimators of E(U) and E(D) are 
not the simple sample arithmetic means.In particular, this 
me true for, 

(1).The log-normal distribution with unknown (to be 
estimated) parameters; often this is used as a nodel 
MOuNdOwn Cl Gepalc tames. In particular, the log 
normal distribution has enjoyed considerable popu- 
larity for representing electronic system repair 
times.Se2 Kline and Almog{[ Ref. 3]. 

(2).the general gamma distribution with unknown shape 
and scale parameters, and 


(3) .the general Weibull distribution with unknown shape 


and scale parameters. 
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B. PURPOSE AND APPROACH 
1. Qbijectives 

Availability estimation by use of the jackknife has 
already been studied by Gaver and Chu.This thesis presents 
two procedures alternative to the above. These nay be 
called the Log-normal Likelihood jackknife (LNLJ) procedure 
and the Moments procedure (MP). In brief, the jackknife 
method has the capacity to reduce the statistical bias of 
estimates of such quantities as system availability, and 
also, and more importantly, to furnish confidence limits 
that behave in a satisfactory manner (economically enclose 
the true availability) despite the fact that underlying dis- 
tributions are unknown. See Millerf{Ref. 4] for a review of 
much of the literature on jackknifing. 

The present estimation procedure makes use of the 
specific assumption of a log normal nodel for repair times 
to compute the maximum likelihood estimate (mle) of avail- 
ability. Properties of the procedure have been studied by 
Monte Carlo simulation. A number of such Simulation results 
are presented in this thesis. Also presented are compari- 
sons with the direct method of Gaver and Chu, and with a 


computationally simple method based on moment estimation. 
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Attention also has been paid to study of the sensitivity of 
MeemuNnLJ procedure to specification error. This refers to 
*he error introduced by assumption of a log-normal distribu- 
ional model and the corresponding maxinun likelihood esti- 
mates when, in fact, another distribution governs the opsear- 
vations. I¢ is found that LNLJ method, which is 
theoretically nost efficient in large samples if the model 
(likelihood) is correct, can sometimes produce consistently 
biased results when the log normal model does not apply. 
Since it 1s nearly impossible to assure log normality fron 
small samples of data, this finding suggests «hat caution is 
mae Order. 

WeeeeLOog nommal Likelihoog Procedure For a Single Unit 

(a) .Data is given as follows, 


U, Un rere cccccvccos ( U (up ines) 


i 


dad, rece ceeceeeee dn (down times) 


xX, =lnd, ,x,=lnd,,......---,x,=lnd, 


(>) .ASSume down times are log-normal.Then, 


E{ D j= exp(u, + lf2e o 47) (cg A 
where 

Bas olen } (Cteoe2) 
ie aa eenD | (1.2.3) 


1 22 





Furthermore, the maxinum likelihood estinates 


mle) of wu, and o? are 
(mie) Wg 5 


ro 3 ca = [ 
a =1/nz in(a.) (1.2.4) 
G2 Vn 8 RR (1.2.5) 
4 i g @ e 


so the mle of E[D]} is, by invariance, 

Sand! ow ~a 

EC DJ = exp(H,+ 1/20 o 2) (i220) 
it is essentially the latter expression that 
will be used in the estimate for availability to 
replace the simple first moment estinate £(D) =d. 


Ie eeieret oo 


5 has been replaced by its unbiased 


version (n replaced by n-1). 

(c).If the above holds, it is advantageous to trans- 
form first (See Mosteller and Tukey{ Ref. 5]j) che 
estimated availability 


3 (0) 


—CeKC_eo_o—X—X—_———— 


B(U) +E (D) 
Pieomueg—lOogGistlec eranstornm of availability is 
In(A / 1-3) = 1nB(u) - (7, +1720 5 2) (1.2.7) 
Jackknifing will be carried out using the sta- 
SUS C 


Lz in (A / 1-R) =1n (T)- (¥+1/2¢5, 2) (in Dees 


ae 





(d) .Recompute a repeatedly, leaving out succes- 


Sively the sample pairs Doe a do ia 


(uid ) 


el) eooeceve y 
Z 


1g 
a5 1/n- om Sy yas (isa 9) 
eset 
x pve ye x. ct Ee Abs Baio 
Sy 200 (x, = x) 74 Ye 2) ate: Wed) 
1a) rayee = 
= - . @ 2 
Lin, 3 in (u) (x 5+ ges Be) (ener) 
Eg c oa meg fee |e ee rr 
(ec) .Compute the pseudovalues as follows; 
PLN, SS, ig ON oe Misip cop ictetety © (lene et 3) 
Recall that L  =L is the result of the con- 
LN LN,all 


puting the quantity to be jackknifed, leaving 


out nene of data. 


(f) .Compute the mean and variance of the pseudova- 


lues, 
- a) 
— D 
PUN W780 0, Bae 
Je! ia 
he | - 2 
a /Nn-1e 2s (P 2 soy, 


(leas 1a) 


(Te2e545) 


(gj) .The jackknifed point estimate of =he availabil- 


ty 1S now; 


A = Spt? ee / leexp (P,, ah 


(h) .Symmetric two sided confidence limits 


dence level (1-21% are derived as 


SGN Soh as mod) eon (i /nes © 
ae See (n-1) SQRT(1/ i ) 
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Follows; 


(eed) 





= B — ® + 
cae Bae t pea (iam )meoon . (177 Sel (ie 2. 18) 


where fy eg 2 (87) is the (1-2 4 /2)100% quantile of 


Student's-t with n-1 degrees of freedom. Then, 


exp (CL) 2xp (CU,, ) 
eer oe AS (Wer 2ien tl 2) 
1+ exp(ch, ) 1 + exp (CU) 


With confidence approximately (1-2 )100%. Note 
that the confidence limits are nearly symmetric 
anounmade ins) uo I/EID]), and not around A. 

(i) .One sided confidence limits at confidence level 
(l-a )100% are derived 
oreo ce (n= 1) SQRT(1/neS_ 2) (1220) 


ef-— 6, * *,_, (n-1) SQRT (1/neS, 2) (ane 1s 


So, one-sided upper confidence limit is 


SS (ieee 2) 
1 + exp (CU) 


and lower confidence limit is 


exp(CL..) 





Ole 22-23) 
1 + exp (CL) 


Note that if up times are exponentially distributed 
then u is actually the mle of ELU], so under the nodel 


assumptions the (transformed) mle of availability is being 
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foreeKkKnif 
meer. 6 ] 
neverthe 
jackknif 


for inve 


ed; this procedure has been validated by Reeds 

for large sampl2s. Simulation studies are 
lass essential for investigating poroperties of the 
2 or alternative techniques for small samples, and 


stigating sensitivity to nodel assumption. 


(a) .Data given as follows 
WU pUspeee eee DL (up times) 
dd, rece ee ed {down times) 
y, =Inu -y, =lnu,,.---,y =lnu 
x =Ind, ,x, =lnd, ,....,x%, =lad, 


(bo) .ASSuMme up and down times are Log-noramai.Then, 


E{ U j= exp (u + WAZeC..<) (AGS 6 1) 
ECD J= exp (¥5+ 1/299 ,2) (ales 32) 
wnere, 

Lys EC lind] ie So 3) 
a ELD | (Aces) 
o*= Gors(eend) (4355) 
fe Var( ind] (453.6) 


Pigemarnore, the mie of u , a4, ua! and a are 
u UW 


_ 
We Sr y Ys (1.807) 
7 jai 
ples ee 
ae 1/ne d. (y, y)2 (eto. 3) 
= nm 
we x =1/ne y F (ie 359) 
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(Cc) 


(d) 


(2) 


(f) 


; | 
o2 = 1 eye 1310 
- /ne 2 (x, ~ x) ( ) 


so the mle of EL U] and ELD] are, by invariance, 


ELUJ= exp(u,+ 1/2¢% 2) ames 11) 
E{ D j= exp (i+ 1/2e 3,?) (1.3.12) 


Pema vOVemnoLas, tne Logistic transform of avail- 
ability is, 

ln (A / (1-A)) = (yt 172 2) - (T+172 62) (1.3.13) 
estimation will be carried ou+ using the statis- 
uC 

P= Tn (A / (1-A) ) =(y+1/2¢s 2) - (x+1/205,2) (1.3.14) 


Simply the point estimate of availability is 


NOW, 

A= exp(P) / 1+exp(P) (Veo. 2)) 
Variance of the statistic P can be calculated as 
FOLLOWS, 


Var (P)=(SE)2= 1/nes,2 + 1/4eVar[s,2) + 1/nes,2 + 


Sg 
\74eVar(s, 2] (1d 600) 
where, 
2%.= - e\2 (1235-17 
Var{s.. j= 1/n C84, (S., )2] (lees ) 
Var[s,? }= T/n Vie (Sc) — | (lee 4 3) 


Syumecanc two sided confidence linits a* conf:i- 
dence level (1-a2)% are derived as follows, 


Pa BO By ee (1.3.19) 


A 





U= p + Z 4/2 (SE)? (ns 20) 
where Z1-a/2t5 the (1- 2/2)100% quantile of the 
Standard normal distribution. Then, 

{exp(L) / itexp(L)} S A S$ {exp(U) / 1+exp (T)} 


(Versa 1) 
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A simulation procedure has been used to compare the 
operating characteristics of LNLJ methodology introducea 
here with other approaches. Specificaily, simulation has 
been used +0 conpute 
(a) the actual coverage of the true availability figure, A, 

by the confidence intervals given by the procedures 

under study, when the nominal coveraga is (1-24)100%, 

(b) measures of confidence interval expected width and vari- 
ance width, 

(c) estimate of the expected point availability estimated by 
the procedure under study. 

The Simulation progran ae welocem, Ln SORERAN LV, and 
the Simulation have been carried out on the I3BM 3033 at the 
Naval Postgraduate School. A mora detailed dascription of 
the program and its major subroutines is given in Apvendix 
A. 


An outline of the simulation procedure now follows; 


% 


Cleese 80 teveed Sample of Up-ctime durations, (u,, i=1, 


fer... 2) amd@e down-tinze durations, (d; t= lg 22ers. en) Were 
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times. The Naval Postgraduate School LLRANDOM package was 
used, along with the International Mathematical and Statis- 
tical Library (IMSL) random number generator. For all cases 
considered (exceptions noted), expected up and down times 
were Ef 0jJ=73.036 and ELD J=3.844231 for the population san- 
pled, so long-run availability A=0.95. Sample sizes of n=15 
and n=25 were simulated, and provided the bases for point 
estimators and for confidence intervals. 

(B) Here are the different distributional situations 
that were sampled. A total of 1000 replications was used to 
evaluate each procedure in each distributional situation. 

(1) Up-times come independently from the iid Expon(,) 
distributions. Down-times come independently fron 
the log-normal(u,c4) distribution. 

(2) Up-times are iid Expon(,). Down times are iid 
Pepe (2) < 

(3) Up-times are iid Expon(i’). Down-times are iid 
Gamma (9,k), being a scale and k a shape parameter, 
so arranged that Ef DJ= Bq? and Var[D]= ({k» y72. A 
Gamma with k>1, 2.g. K=2, represents data tha*+ is 
more tightly grouped around its nean than is true of 


exponentially distributed data; it roughly resembles 





(4) 


(3) 


(6) 


data that is log-normal. A gamma with k<1, e.g. 
k=1/2, represents data of relatively extreme disper- 
Sion and positive skewness as compared +9 exponen- 
tial data having a very long right tail, compensated 
for by a high density near zero. A gamma distribu- 
fronmwath imeeger k (k=1,2,3,....) is often called 
Erlang; realizations are easily simulated by summing 
k independent Exponential realizations. 
Up-times are iid Expon(A). Down-times are iid 
long-tailed h distribution: 

d=E[ D Jeweexp(hew)e(1-h)2 , O<h<1, 
where w comes from Expon(1). This long-tailed h 
distribution possesses no closed-form representa- 
tion. However, it can be shown *o have characteris- 
tics similar to those of the Exponential for small 
values of h, but to have a systematically longer 
EBrgucetall than the exponential. 
Up times are iid long-tailed h distribution. Down 
times are iid Expon (nu). 
Up times are iid Expon(A). Down times are iid Wei- 
pi Ld, Kk), being a scale parameter and k a shape 


parameter, so arranged «that SCD j=91/Ke Pa = Fe te an 


on 





(7) 


Var{D]=o-/“ef r (2/k + 1)- p2(1/k + 1) J. A Weibull 
With k>1, e.g. k=2, represents data that is more 
Eagieny GLewuped shound 1ts mean than 1s true of 
exponentially distributed data; it quantitatively 
resembles data that is log-normal.aA Weibull with 

kK<{ 1, @.g. k=1/2, represents data of relatively 
extreme dispersion and positive skewness as ccmpared 
to exponential data a very long right tail, compen- 
sated for by a high density near zero. 

Up times are iid Expon(A). Down times are iid 
long-tailed log-normal h distribution. See Appendix 


B. 





El... Wars. S 


The methods described were simulated for various @istri- 
butional assumptions which were defined in the Chapter 2. 
Simulation results for each method and set of distributional 
assumptions are shown in Table 7 and Table 2. 

In general, MP has very high coverage factor and satis- 
factory point availability estimation. However, both average 
and variance of confidence length are higher than JK and 
Maboeprocedures.In addition, no consistent superiority of 
LNLJ procedure over JK procedure has been noted. 

Especially under the Exponential up and Gamna (k=1/2) 
down, and Exponential up and Weibull (k=1/2) down times 
assumptions, a consistently low coverage factor and peint 
availability estimation, and a consistently high average 
confidence length and variation have been obtained from LNLJ 
procedure. In these cases, distributional characteristics 
memeaownm times after log-transform provide highly left skewed 
Shape and large variance, and, also these features reflect 


in «he pseudovalues which, in turn, causes left sk 


‘D 
x 
© 
Qu 
W 
a 
rey 
Ae. 
© 


rather than a symmetric shape for the pseudoviius 


Seestribution. 


Zo 





In order t9 compensate for this situation, the Biweight 
procedure, see{Ref. 7], and the Winsorizing procedure, see 
(Ref. 8], were implemented directly on the vseudovalues. 
Results are shown in Tables 3, 4, 5, 6, 7, 8. These proce- 
dures are robust/ resistant methods for establishing confi- 
dence limits on means; they tend to down-weight extreme 
observations that appears as outliers. 

For small sample size, (i.e. n=15) Winsorizing at g=2 
and g=3 levels provide an adequate coverage. But, the poin 
availability estimate and average confidence length remain 
inadequate. In general, the Blweight statistical procedure 
caused reduction of the coverage factor, and it is not yet 
considered effective for compensation. 

Next, a grouped jackknife procedure was applied to these 
cases.Simply, up and down times data were grouped two by 
two, and average of these groups have been taken. Obviously, 
this process reduced the sample sizes to one half and caused 
loss of degrees of freedom. However, smaller variation and 
more stable results were anticipated.Simulation results are 
presented in Tables 9 and 10. Very gocd coverage factors 
and point availability estimates were obtained, but, average 


confidence length remained at a rather high levei. See 
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Meeendix C for histograms of both cases for up and down 
times, and log-transform of down times and pseudovalues 
fter grouped jackknife procedure appiied. 

Other critical simulation results were obtained from 
exponential up and long-tailed log-normal h down +imes 
case.For both JK and LNLJ procedures, above statistical 
procedures were applied, and results were shewn in Tables 
11, 12, 13, 14, 15, 16.In either case, after implementation 
of these statistical representations, the coverage factor 
decreased. Also, see Appendix C for histograms of this 


case. 
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Table 1; Simulation Results for Distributional Assumptions 
Sample size=15 

Underlying Average Variance Point 

Distribution Coverage Width Width Availability 
A. Exponential up JK Q.9380 Ves sia 0. 00 22 ORS 
Log-normal down LN 0.9443 Qewese | 0. 0016 0.9459 
Us! (ls Sia) hte} Og 0. 00 26 Orr S56 
B. Exponential up JK 0.9470 Deve50) 0.0013 0.9471 
Exponential downLN 0.9417 Ore 7Ss. 020129 Oo els A 
MP 0.9868 Oem oom 0. 0196 0.9421 
C. Exponential up JK 0.9423 0.0706 0.0007 0.9465 
Gamma down LN 0.9419 0.0764 0.0010 0.9437 
k=2 yle) NO Bhe! Sil0, 0.1001 0.0031 9) EMSIS Bie 
D. Ditto JK 0.9450 Otis 1 30.0039 0.9442 
k=1/2 LN 0.9161 0.4574 0.0932 0.8064 
Ge 0. 2873 Ome 200970 Ones 
Beeeaxponential up JK 0.9256 Oe io. 0097 0.9437 
Long-tailed h LN 0.9483 eto Or O20 197 O92 87 
h=0.2 MP 0.9780 0.1951 0.0244 0.9412 
EF. Daxct oO TK a0 to 3 Oe tsios.” 0 0343 029357 
h=0.4 LN 0.9406 OO, 2184) 0.0242 0.9374 
MO. 25 13 Op 20S 202-0324 Poe424 


26 





Ge 


H. 


Es 


J. 


K. 


Long-tailed h upJK 


Exponential downLNn 


h=0.2 


D2tto 


Exponential up 
Weibull down 


k=2 


DEtto 


k=1/2 


Exponential up 
Long-tailed log 


normal h down 


he 


JK 


LN 


we 


JK 


LN 


oe 


JK 


LN 


es 


JK 


MP 


JK 


LN 


We 


Ome 55 


0.9058 


V2 7a03 


OPS ic) SW 


0.8552 


0.9437 


0.9476 


0.9449 


Dis ghey Bie 


0.9742 


0.9336 


0.8968 


DG SA As 


Gag1 03 


0.8794 


022009 


0.9001 


27 


Como? 7 


Gras 0 


Gee 73 


0.1470 


0.2100 


0.2377 


0.0641 


0.0672 


O02 17 


0.1764 


UGS) 1 a4 


0.5819 


0.1385 


Criiz 7 


0.117113 


051978 


0.1394 


Cee 


0.0024 


0.0150 


OP 


0.0045 


OPE Ae, 


0.0244 


0.0005 


0.0007 


0.0023 


OS es Be, 


0.0921 


0.0966 


0.0220 


0.0091 


0.0058 


0.93505 


Bg a) 1 site 


Or 0 Od? 


O239357 


0.9163 


0.9284 


0. 9483 


0.9463 


C2950 


0.9404 


0.7649 


0.8464 


0.9416 


0.9483 


OS Isis 


0.9394 


0.9483 


Us 955 





Table 2: Simulation Results for Distributional Assumptions 
Sample size=25 

Underlying Average Variance Point 

Bestribution Coverage Width Width Availability 
Beexponential up JK 0.9436 0.0594 0.0007 0.9470 
Log~normal down LN 0.9473 O70 a2 0. C004 0.9471 
MP 0.9740 Of. 0713 70.0008 0.9571 
B. Exponential up JK 0.9473 020597 0.0004 0.9u71 
Exponential downLN 0.9228 Or0o7 85 0.0052 Om gs22 
MP 0.9945 Oe 202s 0.0086 0.9450 
Meexponential up JK 0.948475 0.0506 0.0002 0.9474 
Gamma down LN 0.9414 Qos 47) 9) 0. 0003 0.9447 
k=2 MP 0.9895 O05 1 0. 0010 Uo Ses 
D. Ditto JK 0.9464 020766 "90.0009 0.9463 
k=1/2 Oe ola 0.4008 0.0790 0.8168 
MP 0.9884 0.4506 0.0841 0.8435 
feeeezponential up JK 90.9313 0.0834 0.0041 0.9456 
Long-tailed h LN 0.9464 Ort 13 07. 0:077 Ors 
h=0.2 MP 0.9867 Oragzo » 0.0108 Omg >0 
6 Dae pelt JK 0.8784 Om2649 0. 0212 0.9421 
h=0.4 IN 02.94 75 One 46 Game 0. 0 132 0.9366 
Wee). 95 93 Oe VS 0am Ois Oils 9 0.9481 
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Table 3: Implementation of "Winso 


and Gamma (k=1/2) 


Coverage 
G=0 0.9103 
G=1 0.9095 
G=2 Oro 27 2 
=3 0.9400 


Table 4: Implementation of "Winsorizing" on Exponential Up 
and Gamma (k=1/2) 


Coverage 
G=0 078256 
G=2 0.38010 
G=3 0.55.36 | 
G=4 One717 


So ed *] 


ee es 


Sample size=15 


Average 


Width 


0.4212 


0.3624 


0.3467 


0.3427 


Varlance 


Wee h 


0.0809 


Ore O5 7-2 


DPMS yeh 7 


ORR0609 


Sample size=25 


Average 


Width 


0.3691 


0.2565 


0.2456 


De SNS) 
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Variance 


Width 


0.0670 


UE isis: 


0.0247 


0.0247 
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Taple 5; Implementation of "Winseorizing™ on Exponentiai Up 
and Weibull (k=1/2) Down Times for LNLJ Procedure 


Sample size=15 


Average Variance Point 
Coverage Width Rialgelse al Availability 
G=0 0179336 U7 22 0.0921 0.7649 
G=1 0.9418 0.5056 0. Ogee Gaas2s 
G=2 Or 95.2.6 0.4699 O20 65)5 Oe siey ay 
G=3 0.9510 0.4464 0.0907 On58'53 6 
Taple 6: Implementation of "Winsorizing" on Exponential Up 
and Weibull (k=1/2) Down Times for LNLJ Procedure 
Sample size=25 
Average Variance Point 
Coverage Width aoe th Availability 
G=0 0. 8456 0.5116 0.0796 027765 
G=2 UeecGs 2 0.3874 0.0516 Os Se7 
G=3 0.9051 His 6 32 O05 13 0.8689 
G=4 0.9316 0.3431 0.0511 0.8829 
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iaole 7: implementation of "Biwelghting" on Exponential Up 
and Gamma (k=1/2) Down Times for LNLJ Procedure 
Sample Average Variance POUR: 
size Coverage width width Availability 
@raginal 5) Ono 03 0.4212 0.0809 Oee2cs 
Biweighted 15 Oe S2i5 05 262.) 0.0448 OSS be lin 7 
Original 25 ORs) ZENS 0.3691 0.0670 0.723.419 
Biweighted 25 0.8786 On lwo 0. 0159 0.9495 


moeete 3: implementation of “Biweighting"™ on Exponential Up 
and Weibull (k=1/2) Down Times for LNLJ Procedure 
Sample Average Variance Point 
slze Coverage width width Availability 
Original 15 0.9336 Oma 5 2 Oe09241 0.7649 
Biweighted 15 0.38423 023153 0.0784 0.9071 
Original ZS 0.3456 02> iS C7076 0. 7765 
Biweighted 25 0.8374 0.2134 0, 6362 02 92s0 





Table 9: Implementation of "Grouped Jackknife" Procedure for 
Exponential Op and Weibull (k=1/2) Down Times 
Sample Average Variance Point 
Method Size Coverage Weed = Width Availability 
JK 14 022209 Dea 136 OE 18 Si ays: 0.9381 
LN 14 0.9486 0.3648 Oe or 7a On Os 
MP 14 0.9390 Vao0 33 0.0672 029073 
JK 24 Ono223 Oewliss9 O20 113 0.9429 
LN 24 059500 0.2469 0.0400 0.9148 
MP 24 Qecos 5 O16 2S is OnOu 52 0.9183 
Table 10: Implementation of "Grouped Jackknife"™ procedures 
for Exponential Up and Gamma (k=1/2) Down Times 
Sample Average Variance Posant 
Method Size Coverage leyobre tek Width Availability 
JK 14 0.9426 OFF tS 310 020030 0.9441 
LN 14 0.9534 E22 66 oO, 036" 0.9243 
MP 14 0.9608 0.2186 0.0419 0.9262 
JK 24 0.9472 07 0839 ORO 04 5 0.9461 
LN 24 Oo 62 0.1448 0.0174 059293 
ie 24 0.9783 Ont? 10 DeO259 On9322 





maple Vi; Implementation of "Winsorizing" on Exponential Up 
and Long-Tailed Log-Normal (h=9.2) Down Times 


Sample size2=15 


Average? Variance Point 
Coverage Width Aidth Availability 

JK 

G=0 0.8968 Oa 1e85 0.0220 0.9416 
G= 1 0.8190 0, 006 Je 0017 0.9544 
G=2 0.7450 07 0585 0.0011 029.516 .0 
G=3 057160 0.054 Orie 209 HESS oS 
LN 

G=0 0.9119 Onaizal 0.0091 0.9483 
G=1 0.8560 0.0898 0.0040 0.9549 
G=2 0.8040 Os 05° 020029 URS Ss 
G=3 0.7460 0.0668 0.0027 On oee 


34 





Meoe 12: Lmplementation of "Hinsorizing" on Exponential Up 
and Long-Tailed Log-Normal (h=0.2) Down Times 


Sample size=25 


Average Variance Point 
Coverage Awdth Width Availability 

JK 

G=0 0.9010 Oze03 2 OeOn2 3 0.9440 
G=2 en a8 0 0.0460 0.0004 0.9554 
G=3 0.6950 0.0416 5 (O00) 0.9563 
G=4 0.6620 Oe0si8 7 0.0003 029568 
LN 

G=0 0.9141 Odie 7 O00 26 On 22107 
G=2 0.8240 0.0594 0.0010 Os 25538 
G=3 Ora SO Ome 3 9 0.0008 0.9573 
G=4 0.512930 0.04388 0.0006 0.9588 
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Table 13: Implementation of "Winsorizing" on Exponential Up 
and Long-Tailed Log-Normal (h=0.4) Down Times 


Sample size=15 


Average Variance Point 
Coverage Azad =n width Availability 

JK 

G=0 Ooo 0.1578 Os Oe SIOiSs Qn J594 
G=1 0.7560 020720 O20 C26 0.9542 
G=2 0.6870 Oe Ol Syeye: OTOUrS OR 010 
G=3 0.6460 020587 OV0012 Or 931615 
LN 

G=0 O22 009 OF ssa 020136 0.9483 
G=1 0.8470 0.1089 0.0069 0.9543 
G=2 0.7800 0.0866 Oh lO) By | Ooo 7 6 
=3 0.7180 020762 0.0048 0.9603 
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Table 14: Implementation of "Winsorizing" on Exponential Up 
and Long-Tailed Log-Normal (h=0.4) Down Times 


Sample size=25 


Average Variance Point 
Coverage Wemce con Width Availability 

JK 

G=0 OmecsS 0.1184 Or Omions 0.9426 
G=2 OmiGid 110 0.0473 0.0007 Omes6 0 
=3 Ome 2.0 0.0411 0.0004 Ong57 2 
G=4 Ono 5 0 Gr.0 350 0.0004 Oegoa7¢ 
LN 

G=0 R003 0.0918 O00, eS yy 
G=2 0.8000 0.0666 0.0018 06 vieitexs 
G=3 0.7660 OPO 565 0.0014 06 SIE SKe) 
G=4 Oo. 700C 020331 0.0012 OG shell) y 
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jmmole 15: Emplementation of "Biweighting" on Exponential Up 
and Long-Tailed Log-Normal (h=0.2) Down Tines 
Sample Average Variance Poent 
size Coverage width width Availability 
JK 
Original 15 OG Soke: O. 12:8 0.0220 0.9416 
Biweighted 15 0.6870 OR ONS) 26 0.0010 0.9554 
Original 25 Dee0n.0 0. 102 Oh p48 0.9440 
Biweighted 25 0265160 8) 10 SNS 7 0.0002 059569 
LN 
Original 15 ee is Omi 7 0.0098 0.9483 
Biweighted 15 OR 710) 16) 0.0604 0.0049 0.9607 
Original 25 0.9741 0.0817 0.0026 029507 
Biweighted 25 Omori 0 0.0406 0.0004 0.9627 
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Gable 16: Implementation of "Biweighting" on Exponential Up 


and Long-Tailed Log-Normal (h=0.4) Down Times 
Sample Average Variance Point 
size Coverage width width Availability 
JK 
Original 15 OQ. S794 Oels 78 0.03105 0.9394 
Bliweighted 15 G26 2 10 0.0499 0.0010 029573 
Ofiginal BS Cesc. Oni ai8 Dae iyi s: 0.9426 
Biweighted 20 Ve 5100 0.0340 OOO Key O92 90 
LN 
Original 15 0.9009 0.1394 0.0135 v2 9483 
Biweighted 15 © lS) SS) 0.0614 0.0025 0.9642 
®eiginal as 0.9003 0.0918 C007 0.9511 
Biweighted 2 O.a264 9.0394 G2 0005 Gageos 
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IV. CONCLUSIONS 


The LNLJ procedure works well when the actual down tines 
are Exponential or Log-normal. However, the orocedure is 
very sensitive to the sample variance. If th2 sample vari- 
ance tends to vary greatly, sometimes becoming 2xcessively 
large, then the method tends +o fail; the intervals shifted 
down. For instance, long tailed Samma, Weibull and Log-nor- 


Mal distributed down times provide extremely large 


" 


ample 


variances, and under these circumtances the method 


ih 
}-46 


ails. 

mPemeciccher study, one might apply jackknifing leavin 
Out one observation at a time, instead leaving out a pair. 
Obviously, this is going tds increase the degrees of freedon, 
but might by help reducing the variance and the confidence 
interval length. On the other hand, grouping more than two, 
Might provide more stable results. 

In addition, in order to see what happens in one-sided 
availability estimation, Exponential up and Loq-normal down 
and Exponential up Weibull (k=1/2) down tines were simulate 


and results are shown in Appendix D. 


40 





APPENDIX A 


Simulation program consists of main program and several 
subroutines. Main program, for JK, LNLJ and MP procedures 
computes availability confidence limits and point availabil- 
ity and, scores the coverage for each replication. Then, 
after 1000 replications computes the statistics of «hese 
parameters and prints out the results for given underlying 


Seestributions. 


|< 


| +e 


Subroutine CONP computes two-sided confidanca limits of 
a given data vector. Subroutine MOMENT computes the fourth 
moment of given data vector. Subroutine LNLJ generates the 
Meena down times irom given undeclying distributions, and 
takes the log transform of these, and computes the means and 
the variances of these data vectors. 

Subroutine BIWGT computes the Biweigh*t estinates of a 
Given data vector, and uses Subroutine MEDIAN for computin 


the median of this data. Subroutine WINSOR Winsorize a 


given data vector for given "g" lavel. 
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(1) Let W be Llog-normal random variable which 
InWaA, N(u,c?)- k moment of W as follows, 
ec a* j= exp (ku + 1/f2ek2 o2), 
see {Ref. 9]. 
SO, 
EC Wj= exp( ut 1/20 2) 
E{ #2 j= exp (2u + 2 9?) 
(2) Create a stretched log-normal random variable 
D= ceWe( i+thewW) 


where c and h are constant and 0<h<1. 


ECD J= c{E{ ¥]+theE[ w? }} 
E(DJ= cfexp(u+1/2e 5g 2) + heexp (2p +2 o2)} 
= ELD) { exp(u+l/2e 52) +theexp(2 p +2 52) }7t 3 (a.1) 
(3) On the other and, we can write the median for which 
is the function of the expected value as follows, 


M= d(0.5)= E[D] 


M= ceexp(uye(Ttheexp( u))= PEL D] (a. 2) 


a2 





here yeascetne Median location factor of the new dis- 
wemouttOn. Lf we put ecuation (a.1) in equation 
fae 2), we obtain 


exp (u) e(1+heexp (u)) 


VeE[D]= E[D] 





exp(u+l/ze 52) + heexp(2y+2 o@) 
After the cancellations, we obtain 
vy {exp(u + 1/2¢, 7) +theexp(2y + 2 o%) J= 1+heexp (u) 
beexp (1/2 a7) - 1 = heexp(y)e(1 -y eexp(2 o?)) 


1- UV eexp(1/2e , 2) 


BI 
i 





exp (-1y) 
ee C2 ig 6) = | 


If h=0, then yp = exp (-1/2e , 2) 


1 + SXp (1p) 





tf h=1, then vy 
Sxapee won sewn (172°) oH) 

So, if we get 

Eee lyceig =) = YS exp{ tl +2 ,*)+exp(i/2e .*) 


then we obtain h between 0 and 1. 


a] 
tw 





The LNLJ method fails for Exponential up and Gamma 
(k=1/2) down, Exponential up and Weibull (k=1/2) down, and 
Exponential up and Long-tailed Log- normal h down times 
cases. In order to demonsrate the distributional charecter- 
istics of this cases, histograms are obtained from sinula- 
[eomecOr up times, down times, log transformations of down 
times and pseudovalues. Also, histograms are repeated after 
grouping the observations in order to observe the grouping 
Semess On the parameters. Distributional charecteristics of 
these parameters such as means, variances, moments and 


shapes can be observed fron these histograms. 
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of LNLJ procedure for Exponential up and Weibull (k = 1/2) down times. 
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QNE-SIDED AVALLABILITY ESTIMATION 
Best and worst cases for LNLJ procedure were simulated 
for one-sided availability estimation, and results are shown 


in Tables 17 and 18 as follows; 


Table 17: Exponential Up and Log-normal Down Times 
N=15 N=25 
Js LN JK LN 
Low. Availability Bound Cover Cee Jeon Oo. 999° 1,000 
Upp. Availability Bound Cover O2333900.099 0.966 0.965 
Avg. Lower Availability Bound Cee Immosrs7 9 6.9190 027908 
Avg. Upper Availability Bound 0.966 0.964 0.966 0.966 
Table 18: Exponential Up and Weibull (k=1/2) Down Times 
N=15 N=25 
JK LN JK LN 
Low. Availability Bound Cover eee ong 90 O29 70.997 
Upp. Availability S3o0und Cover HaocomOncos 9.955 0.705 
AVg. Lower Availability Bound Ons5ue eo. 460 9.9381. 0.500 
Avg. Upper Availability 3o0und Geom oes 0 0.975 0.955 


vy 





Notice the low Average Lower Availability Bound provided 
by the LNLJ procedure for the Weibull (k=1/2). S8esults show 
that the one-sided availability estimation behaves as the 


two-sided availability estimation. 


78 





Deo ReGen er ERENCES 


Swedish Board for Been ence) DG Goons STU-Rep 
e yy 


75-4090, Test >? y J. 
and QO. Sjorktun , P.wt-S, bDeCcemasr $78. 


Given Dee r., ond Chu, 3. 58., “Jackknife Estimates 
of Component and System Availability," at 
trics, YVol.z1, No. 4, p. 443-449, November 9°79, 
ao eos and cuted, 2., Sujeability of Th 
LO Deo e aoc On fOr vepas *anes.  p 
Sented at {hternacional Conrerence on rReiia 
ane NWaeeueate negli ty, Lannion, rrance, 8-12 


—~normal 
















Miller, R. G., “rhe Jackknife - a 
PWaeerweevol. Of, HO. 1, p. 1-15, 7 


Mosteller, F., and Dakey Jed. Bate Analysic. — 
Including Sn bey es fapter Lm Manclbook.o Cia 
Sycology, Vol. 2, aAdd@ison- Wesley , 1968- 7 

Reeds, T., "On Jackknifing Maximum Lixel: 


— 


: 
feweOr UeNnna oeOr StaesSstics, No. 6, Pp. 7 


1978 

Mosteller, F., and Tukev, J. W., Data Analvsis and 
Hee GEE p. 205, 358-354, Addzson- Nes Sy , 
Dixon, W. 
tatistical 






wommson, Wee b., and Kotz §S Continuous Univariate 


¢ . * @ 
mee Vol. 1, p. 115, John Wiley and sons 
Rich, : 


iS, 





10. 


eves 


Cae 


2 


ee 


PEA VESMRLBUTICN LIst 


Defense Technical Information Center 


Gamneror Station 
Alexandria, Vircinia 22314 


Library, Code O1&%2 
Naval Bost ae School 
Monterey, fornia 93940 


Department Chainman, Ccoce 35 
Department of Operations Research 
Naval Postgraduate School 
Monterey, California 93940 


Professor D. P. Gaver, Code 55GV 
Department of Operations Research 
Naval Postgraduate School 
Monterey, California 93940 


Bogazici Universitesi 
Yon Eylen Analizi Bolumu 
Besiktas, Istanbul, Turkey 


DZpe reas 2 
Bakanliklar, Ankara, Turkey 


Weomeeemiarpokula K.ligz 
Heybeliada, Istanbul, Turkey 


Ortadogu Teknik Universitesi 
Yon Eylem Analizi Bolunu 
Ankara, Turkey 


Deedee, A. N. Fiaathen 
any 
1 Loddef jord, Norway 


= uated oe bfelieos 
2 Malloway Lane 
Monterey, California 93940 


toms. LaVvuz 
Baqlarbasi Mah. 
emg Topel Cad. 
5. SOk No.1 
Kirikkale, Turkey 


Lt. Davut Kirca 

Ipras Sitesi be Sok No.4 
Mieudes=*iik, Kocaeli, Turkey 
Professor NaetoG sest, Code S5SFO 


Department of Operations Research 


R 
Naval Dostd 
Monterey, ¢ 


LT Barbaros Aba 

Denzel tierilosu K.ligi 
Peale tila Komodor lugtt 
Golcuk-Kocaeli-TURKEY 


80 


No. 


Copies 


NO 

















Thesis 
ALI35 
Cel 





~— pur te g 
Cho awe 


Investigation of 
alternative methods 
including jackknifing 
for estimating point 
availability of a 
System. 


Aba 





